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Abstract 

We constructed linear algorithms of sampling recovery and cnbature formulas on Smolyak grids 
parametrized by m £ N of periodic d-variate functions having Lipschitz-Holder mixed smoothness a > 0 
based on B-spline quasi-interpolation, and studied their optimality. We established lower estimates (for 
a < 2) and upper bounds of the error of the optimal sampling recovery and the optimal integration on 
Smolyak grids, explicit in d, m and the number v of active variables of functions when d and m may be 
large. 

Keywords and Phrases High-dimensional sampling recovery • Linear sampling algorithms • Smolyak grids 
• Lipschitz-Holder spaces of mixed smoothness ■ B-spline quasi-interpolation representations. 

Mathematics Subject Classifications (2000) 41A15 • 41A05 • 41A25 • 41A58 • 41A63. 


1 Introduction 

We are interested in sampling recovery and cnbature of functions on 1-periodic at each variable. It is 
convenient to consider them as functions defined in the d-torus = [0, l]'^ which is defined as the cross product 
of d copies of the interval [0,1] with the identification of the end points. To avoid confusion, we use the notation 
to denote the standard unit d-cube [0,1]'^. 

For m G N, the well known periodic Smolyak grid of points G'^{m) C is defined as 
G‘^{m) := U = 2“'=s :kGN‘^, |fc|i = m, s € 

where I‘^{k) := {s G Z(|. : sj = 0,1,..., 2^^' — 1, j G [d]}. Here and in what follows, we use the notations: 
Z+ := {s G Z : s > 0}; xy := {xiyi, ...,Xdyd)] 2“ := (2^i ,..., 2“'^); |a;|i := Yh=i 1**1 Ur G M denotes the 
set of all natural numbers from 1 to n; Xi denotes the ith coordinate of a; G i.e., x := (a;i, ...,Xd)- 

The sparse grids G‘^{m) for sampling recovery and numerical integration were first considered by Smolyak 
[50] . In [5T]-[53] and [H]-[I0] Smolyak’s construction was developed for studying the sampling recovery for 
periodic Sobolev classes and Nikol’skii classes having mixed smoothness. Recently, the sampling recovery for 
periodic Besov classes having mixed smoothness has been investigated in ES1|55]. For non-periodic functions 
of mixed smoothness linear sampling algorithms on Smolyak grids have been recently studied in |34j (d = 2), 
[HHS], using the mixed tensor product B-splines. Methods of approximation by arbitrary linear combinations 
of translates of the Korobov kernel on Smolyak grids of functions from the Korobov space which 

is a reproducing kernel Hilbert space with the associated kernel have constructed in m- 
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In numerical applications for approximation problems involving a large number of variables, Smolyak grids 
was first considered in |3f)j in parallel algorithms for numerical solving PDEs. Numerical integration on Smolyak 
grids was investigated in |19) . For non-periodic functions of mixed smoothness of integer order, linear sampling 
algorithms on Smolyak grids have been investigated in [2] employing hierarchical Lagrangian polynomials multi¬ 
level basis. There is a very large number of papers on Smolyak grids and their modifications in various problems 
of approximations, sampling recovery and integration with applications in data mining, mathematical finance, 
learning theory, numerical solving of PDF and stochastic PDF, etc. to mention all of them. The reader can 
see the surveys in [21 [251 [20] and the references therein. For recent further developments and results see in 
[231122121 Hill]. 

The Smolyak grids G^(m) and their various modifications (in partucular, the grids G{m) and G''{m)) are 
very sparse. The number of knots of G^{m) is smaller than and is much smaller than 2‘^'", the 

number of knots of corresponding standard full grids. However, for periodic functions having mixed smoothness, 
they give the same error of the sampling recovery on the standard full grids. 

Quasi-interpolation based on scaled B-splines with integer knots, possesses good local and approximation 
properties for smooth functions, see 0 p. 63-65], 0 p. 100-107]. It can be an efficient tool in some high¬ 
dimensional approximation problems, especially in applications ones. Thus, one of the important bases for sparse 
grid high-dimensional approximations having various applications, are the Faber functions (hat functions) which 
are piecewise linear B-splines of second order [2 Uni dl dl [Ml H HI- The representation by Faber basis can 
be obtained by the B-spline quasi-interpolation (see, e. g., HSj). 

The object of our interest in sampling recovery are functions having Lipschitz-Holder mixed smoothness 
a > 0. Let us introduce the space i7^(T‘^) of all such functions. For univariate functions / on R, the rth 
difference operator A]] is defined by 

1=0 

If u is any subset of [d], for multivariate functions on the mixed (r, rt)th difference operator is defined 

by 


A 


r,u 

h 


Hal. Kf = i. 

i^u 


where the univariate operator A]], is applied to the univariate function / by considering / as a function of 
variable Xi with the other variables held fixed. If 0 < a < r, we introduce the semi-norm |/|//a for functions 
/ G G{T^) by 

|/|^.(„) := sup n^““ll^r(/)llc(TQ 


(in particular, (T<i)(0) = ll/llc(T<i))- The Lipschitz-Holder space of mixed smoothness a is defined 

as the set of functions / G G(T'^) for which the norm 


||/||^^(Td) sup 

uC[d] 

is finite. The non-periodic space can be defined in a similar way with a slight modification. 

In what follows, for the space the parameters a and p are fixed and therefore their value is 

considered as a constant. Denote by the unit ball in Paralleling with we will consider the 

sampling recovery of functions from its subsets and , 1 < v < d. The set consists of all functions 
/ G such that f{x) =0 if xj = 0 for some index j G [d]. It can be seen also as the subset of the unit ball in 
the space with zero boundary condition. The set consists of all functions / G having at most iz 

active variables. Here, we say that the jth variable Xj is active for a function / on T”* if for all i ^ j there are x* 
such that the univariate function g{t) = /(a;f, ...,x*_l,t,x*_^_l, ■■■,x’^) is not a constant as a univariate function 
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in t. With this definition we have = U^. The set U^’' can be a model of the objects in a d-variate space 
depending only a few v (much smaller than d) of variables without in general, exact information about them. 

For sampling recovery of functions / € we use the linear sampling algorithm 

on the Smolyak grids 

G{m) ■- {£, = 2-'=s : fc G N^, \k\i = m, s e 

where 4^ = {¥’4}{eG(m) ^ family of functions on := {s G : Sj = — 1, j G [d]}. For 

sampling recovery of functions / G , we use the linear sampling algorithm 

on the Smolyak grids 

G^{m) := = 2“^s : k G N'^, |A:|i = m, | supp(A:)| <v,s& I'^{k)}, 

where is a family of functions on T^, /^(fc) := {s G : sj = 0,1,..., 2^^' — 1, j G [d]}, 

supp(fc) denotes the support of k, i.e., the subset of all j G [d] such that kj ^ 0, and |A| denotes the cardinality 
of the finite set A. 

Let us notice the following. The number of points in the grid G{m) is 

|G(m)| = E 11(2'=^ -1) < (1-1) 

feeN"^: |fc|i=m i=l 2 

and the number G'^ (m) of points in the grid G'^ (m) is 

iG-wi = E E n < 2”(()(7ri').(12) 

\u\<i> supp(fe)=u, |fe|i=m jGsupp{k) 

The choice of the Smolyak grid G''(rn) for sampling recovery of functions / G is quite natural since in 

general, their active variables are unknown. Moreover, the grids G'^(rn) polynomially depending on d, is full if 
and only li m > v, and the grids G(m) are empty if m < d. For this reason, we are interested in constructing 
linear sampling algorithms on the grids G'^(rn) and G{m) for m > v and m> d, respectively. 

In high-dimensional approximation applications, knowing the number of sampled function values on the 
Smolyak grids G''{m) and G{m) as in (11.11) and (jl.2p and linear sampling algorithms 5m(^m, /) and /), 

one attempts to estimate the errors of ||/ — 5m($m, /)||p and ||/ — f)\\p as a function of three variables 

m, d, v when m and d may be very large. To study the optimality of these algorithms let us introduce the 
quantity of optimal sampling recovery Sm{Fd)p on Smolyak grids G{m) with respect to the function class Fd by 

^m{,Fd)p inf sup 11/ f)\\p^ 

and the quantity of optimal sampling recovery s'^{Fd)p on Smolyak grids G'^{m) with respect to the function 
class Fd by 

s';,iFd)p := mf sup ||/- /)||p. 


In traditional estimation of the error, a typical form of lower and upper bounds of Sm(^(S>)p ^ad Sm{U^)p 


IS 


c'(d)2-“-m''-i < Sm{U^)p < st{U^)p < c(d)2 


-am^d-l 


(1.3) 
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where c(d) and c'(d) are a constant not explicitly computed in d (see, e.g., [T31 HZl [33] ). In high-dimensional 
approximation problems using function values information, the parameter m is main in the study of convergence 
rates of the approximation error with respect to m going to infinity. However, the dimension d may hardly 
affect this rate when d is large. In the present paper, we establish upper and lower bounds for Sm{U^)p 
and s'^(U^‘^)p explicitly in m and v,d as a function of three variables m, d, v. Notice for example, that in 
traditional estimates the term is a priori split from constants c((i, a,p) which are actually a function 

of dimension parameter d, and therefore, any high-dimensional estimate based on them may lead to a rough 
bound. By a combinatoric argument, it is more natural and suitable to investigate upper and lower bounds for 
Sm{U^)p and in the form 

(7'(d)2—(1.4) 
C'(z.)2—<s:;(C/“a < (1-5) 

or in the more refined form 

i'(m, d) 2-“™ < Sm{UZ)p < A{m, d) 2““™, (1.6) 

A'im, v) 2-“- < sl,{U^np < A{m, v) 2““-. (1.7) 

In the present paper, we estimate Sm{U^)p and s‘^{U^'^)p in the forms (11.31) . (I1.4I) - (I1.5I) and (ll.6l) - (ll.7L 
stressing in finding lower and upper bounds for the accompanying ’’constants” explicit in d or in ^ or in m, v, 
respectively. To obtain these upper we construct linear sampling algorithms on the Smolyak grids G{m) and 
G''{m) based on a B-spline quasi-interpolation series especially, on Faber series related to the well-known hat 
functions. Due to the complicatedness, we are restricted to give these lower bounds only for the case a < 2. As 
consequences we obtain upper and lower bounds for the quantities of optimal cubatute formula Intm(17^)p and 
Int((j(17^’^)p explicit in m and d, v. 

Related to the problems investigated in the present paper, is the problem of hyperbolic cross approximation 
of functions having mixed smoothness in high-dimensional setting in terms of n-widths and e-dimensions which 
have been investigated in [31111137]. 

The paper is organized as follows. In Section H we establish upper and lower bounds and construct linear 
sampling algorithms on Smolyak grids based on Faber series for Sm{U^)p and s‘^{U^'')p for 0 < a < 2. As 
consequences, we derive upper and lower bounds for the error of cubature formulas on Smolyak grids and of 
optimal integration on Smolyak grids, and lntm{U^)p and lat'^{U^'^)p. In SectionH we extend the results on 
upper bounds of s'^{U^'')p and Int((j(C/^’'^)p obtained in SectionUto the case of arbitrary mixed smoothness 
a > 0, based on high-order B-spine quasi-interpolation representations for function from In SectionH 

we give some example of polynomials inducing quasi-interpolation operators. 


2 Sampling recovery based on Faber series 

2.1 Faber series 

Let M 2 {x) = (1 — |x—l|)+,a;Gl, be the piece-wise linear B-spine with knot at 0,1, 2 (the hat function), where 
x+ max(a;,0) for x S R. Since the support of functions M2(2^+^-) for /c S Z+ is the interval [0, 2“^] C I, 
we can extend these functions to an I-periodic function on the whole R. Denote this periodic extension by pk- 
The univariate hat functions ipk,s are defined by 

(/9o.o := 1, ‘Pk,s ■= ^k{- - 2s), fc > 0 s e Z{k), 
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where Z(0) := {0} and Z{k) := {0,1,..., 2^ ^ — 1}. Put Z'^{k) := ^ ^ ^ Z‘^{k), define 

the d-variate hat functions 

d 

‘Pk,s{.x) ■= W_‘^ki,Si[Xi), 

and the d-variate periodic Faber system J^d by 
Fd-.= Wk,s-s&z\k), fcez^}. 

For functions / on T, we define the univariate linear functionals \k,s by 
Am(/) := -iA^_.(/,2-'=+is), fc>0, andAo.o(/) := /(O). 

Let the d-variate linear functionals Xk,s be defined as 

Afe,s(/) := (Afc2,S2 (■•■Afe,j,Sd (/)))> 

where the univariate functional Xki,si is applied to the univariate function / by considering / as a function of 
variable Xi with the other variables held fixed. 

Lemma 2.1 The d-variate periodic Faber system J-^ is a basis in C{T‘^). Moreover, a function f G C(T^) can 
be represented by the Faber series 

f = Y. ^kif) = Y Y Xk,s{f)Tk,S'> (^-i) 

feGZ{ fceZ{sGZ‘^(fe) 

converging in the norm ofC(T‘^). 


Proof. For the univariate case (d = 1), this lemma can be deduced from its well-known counterpart for non¬ 
periodic functions on I (see, e.g., [551 Theorem 1, Chapter VI]). For the multivariate case (d > 1), it can be 
proven by the tensor product argument. □ 

Put := (fc e : supp(fc) = u} for u C [d], and use the convention a;° = 1 for a; S [0,oo]. 

Theorem 2.1 Let 0 < p < oo and 0 < a < 2. Then a function f G id^(T^) can be represented by the series 
m converging in the norm of C{T‘^). Moreover, we have for every k G Z{.(m), 

||%(/)||p < 2 -I“I(p+1)-I“I/p2-“I'=Il (2.2) 


Proof. The first part of the lemma on representation and convergence is in Lemma l2.II Let us prove ( 12 : 211 . We 
first consider the case p = 00 . Since 0 < J2sGZ’‘{k) Tk,s{x) < 1, we have for every k G Z{.(m) and x G T'^, 


\qk{f){x)\ 


sup 

|Am(/)I 

Y Tk,six) < 

sup 

|Am(/)I 

sGZ^{k) 


sGZ'^(k) 

sGZ^{k) 


sup 

|n[- 


1 

VI 

l“lcc“(/,2-'=) < 2-l“l2-“l'=b|/|^„(„) 

s€Z^{k) 

jeu 

z 
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This proves (1^ for p = oo. Hence, if 1 < p < oo, we have also that for every k € Z'^{u), 


dx 


hk{f)\\p = [ >^k,s{f)‘fk,s{x) dx < sup |Afc,s(/)|P f ^ Pk,s{x) 

seZ‘i(fe) sGZ’i{k) JT’i g^zd-{k) 

= sup \Xk,s{fW\Z‘^{k)\ [ \(pk,o{x)\P dx 

s&Z'^{k) 

2 — ^3 

= sup \XkAf)\^\Z\k)\r[2 [ \2'^^xAdxj 
sez-^ik) Jo 

< \f\Hg^(u)Y [2-l“l2l'=li] [2l“l(p + 1)-I“l2-I'=li] 

= 2-l“l(p+l)-l“l2-^’“l'=lM/|^.(„). 

This proves (1^ for 1 < p < oo. Finally, the case 0 < p < 1 can be proven in a similar way starting from the 
inequality 

hkif)\\^ < \^k,sifW [ \<Pk,s{x)\P dx. 

seZ'^ik) 

The proof is complete. □ 


2.2 Auxiliary lemmas 

For m,n gN with m > n, we introduce the function Fm^n ■ (0,1) —> K. by 

s=0 ^ '' 

From the definition it follows that Fm,n{t) = (^) F and conseqnently, 

< Fm+l,nit), t e (0, 1). 


We will need following equation proven in [21 (3.67), p.29]. 


1 ” 

Fm,nit) = ’YZTf F! 


s =0 


1 - t 


For nonnegative integer n, we define the function 


(l-2t)-i. 

t < 1/2, 

2(n + l), 

t=l/2. 

(2t-l)-i[t/(l-t)]"+i, 

t > 1/2. 


Lemma 2.2 Let m,n G N, m > 2n, and 0 < t < 1. Then we have 


Fm,n{t) < 


m 


bnit). 


Proof. We have 



< (n + 1) 



(2.3) 


(2.4) 


(2.5) 
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Hence, by (12.41) the case t = 1/2 of the inequality (12.51) is proven. Next, we have for a: > 0 and x ^ 1, 




/m\ x"+^ - 1 

\n J X — 1 


and 


x"+i - 1 
X — 1 


(x — 1) ^x"+^, 
(1 -x)-i. 


X > 1, 
X < 1. 


By using the last two inequalities for x = {1 — t)/t from (12.41) we prove the cases t < 1/2 and t > 1/2 of the 
inequality (12.51) . □ 


Lemma 2.3 Let to, n £ N, to > n, and 0 < t < 1. Then we have 
1 ^ t 




1-t \ l-t 


1-r 

exp) —-— I TO" 


Proof. We have for x > 0, 

n / \ n 


m\ 


^ E"' 7! ^ 


s!(m — s)! 

Using this estimate from (1^ we deduce the lemma. □ 


2.3 Upper bounds 


In this subsection, we employ the representation by the Faber series (ED to construct linear sampling algorithms 
•) and •) for functions from 17^ and respectively. By use of Theorem 12.11 and auxiliary 

lemmas in the previous subsection we establish upper bounds for the error of the sampling recovery by these 
algorithms and therefore, for Sm{U^)p and 

We first construct linear sampling algorithms functions from U^. By the definition we can see that 



For TO £ we introduce the operator Rm for / £ JJ^ by 

RM) := i: «(/) =1:1; 

/cGN^:|fc|i<m s^Z^{k) 

The operator defines a linear sampling algorithm /) on the grid G{m) by 

Rmif) = = Y. 

^eG(m) 

where ■= 

d 

C = 2“''s, fc £ N"*, s£l^{k), 
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and the univariate functions ipk,s, k s € are defined by 


'02-''s 


1 — ¥^0,0, 

^0,0, 

< 

' fk,j, 


k = 1, s = 0 
k = 1, s = 1 

fc > 1, s = 2j, 0 < j < 2'=-! - 1, 
fc > 1, s = 2j + 1, 0 < j < 2'=-! - 1. 


For 0 < p < oo, a > 0 and m > I — 1, put 
b = bia,p):= 2(2“-l)(p+l)i/^ 


1 — 1 / \ 

I3{l,m) = f3{a,l,m):= (2“ - 1)“'^ f ^ j (2“ - I)'*, /3(0,m):= 1, 

s=0 / 


and 


1 , a> 1, 

a°{d) = a°{a,p,d):= |2“ - 2|-l[2(p + 1)1/^]-'" x <( d, a = 1, 

a<l. 


The following theorem gives upper bounds for the error ||/ — Rmif)\\p for / € and therefore, for 


Theorem 2.2 Let 0 < p < oo and 0 < a < 2. Then we have for every m > d, 
SmiUZ)p < sup ||/-i?„(/)||p < [2(p+l)l/P]-i/3(d,m)2-“- 


( 2 . 6 ) 


< exp(2“ - 1)6“''2 


— d r\ — am ^d—1 


Moreover, if in addition, m > 2{d— 1), 

°Sm{U^)p < sup \\f-R^{f)\\^ < a°{d)T 


f&us, 


d-l 


(2.7) 


Proof. Let us prove the lemma for 1 < p < oo. It can be proven in a similar way with a slight modification for 
0 < p < 1. From Theorem 12.11 (12.4|) and Lemma o it follows that for every f G and and every m > d, 

\\f-Rm{mp< E 2-'^(p+i)-"/^2-“i'=i^ 

= 2-'^(p+l)-'^/P 2-“l'=li =2-'^(p+l)-'^/P^ 

feeN‘^:|fc|i>m i=l \ ) 

= 2-‘^ (p + l)-'^/PF„,d-i(2"“) = [2(p + l)i/p]-d 2-“’"/3(d, m) 

< exp(2“ - 1) b-^ 2"“™ 

The inequalities (EH) are proven. The inequality (EZD can be derived from (12. 8|) by applying Lemma 12.21 for 
t = 2-“. □ 




Notice that some upper bounds of ||/ — Rm{f)\\p with p = 2,oo for functions with zero boundary condition 
from the Sobolev space were established in [5]. 


We next construct linear sampling algorithms functions from . The definition of implies that 
= {/eC/“ : 3uC[d], \u\ = v. f = Y. Y. 9fc(/) 

For m € Z+, we define the operator R'^ by 

R'LU) ■■= E = E E 

: I supp(fc) I \k\i<m : | supp(A:) | <1^, |fc|i<m sGZ'^{k) 

For functions / on R’^ defines a linear sampling algorithm /) on the Smolyak grid G'^{m) by 

RW) = = E 


where := {t/'^}^eG‘'(m)- Notice that the construction of the operator R'^ is quite reasonable since in our 
setting we assume that the active variables of a / S are unknown. 

Put 


7 ( 1 /, to) = 'y{a,p,iy,m) := E ^ *(P+ 1 ) 

1=0 ^ ^ 


and 


{ 1, a > 1, 

0 = 1 , 

(2“-!)-^', a<l. 

Theorem 2.3 Let 0 < p < 00 , 0 < a < 2 and 1 < ly < d. Then we have for every m > v, 

< sup \\f-RW)\\p < 


< exp( 2 “ - 1 ) ( 1 + !/&)“'2 


y rt—otm ^y — 1 


Moreover, if in addition, m > 2{y — 1), 


< sup \\f-R'L{f)\\p < a{v)2- 


v-l 


(2.9) 


( 2 . 10 ) 


Proof Let us prove the theorem for 1 < p < 00 . It can be proven in a similar way with a slight modification for 
0 < p < 1. Put t = 2““. We first consider the case v = d. From Theorem 12.11 (12.41) and Lemma 12.31 it follows 
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that for every / G and and every m > d, 

ii/-<(/)iip < E E 

tiC[d] /ceZ^(u): |fe|i>m 

< E E 2“l“l (p+ 

uC[d] kGZ^iu): \k\i>m 


0 /cGN^: \k\i>m 

= E (I) ^^‘(p +') 2 -"<”+» 

/^O ^ ^ j=l ^ 

Further, applying Lemma [231 to f|2.1ip we get 

||/-<(/)||p < 2 -“(-+i)E(?) 2 -'(^^+ 1 ) 

1=0 ^ ' 




1 / t 




1 -1 
exp I - I m' 


i-i 


< 2-“(-+i) exp(^l^^ E [fj 2-' (P + 1)”'^^ 

= 2 -“™exp|^ 2 “-l^m^-iE 2 -'(p+l)-'/P |^ 2 “ - 

= exp( 2 “-l)(l + l/ 6 )'^ 2 -“™m‘^-^ 

The inequality (12.911 is proven. 

Let us verify (I2.1()L Applying Lemma [2.21 to Fm,i-i{t) in (12.1111 we have for m > 2{d — 1), 

||/-<(/)||p < 2-“(-+i)E(;)(;ri)2“'(P+ir'^"^*-i(2-“) 

< 2-“(™+i) (rf™ i) E (^) 2”' (p + bi-i(2-n 


( 2 . 11 ) 


This completes the proof of the theorem for the case v = d. 

We now consider the case v < d. Let / G . Without the generality we can assume that the active 
variables of / are among Xi, ...,x^. Hence, we have 

= E E <lk{f){x), R’^if) = E E dk{f){x). 

u(l[iy] kGl,^{u) ttCM fceZj(ii): |fc|i<m 

Considering / as a function in U^{T'') and i?((, as an operator in C(T'^), and applying the proven case d = v 
to ||/ — .Rm(/)||p, we prove Theorem 12.31 □ 
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2.4 Lower bounds 


Theorem 2.4 Let \ < p < oo and 0 < a < 2. Then we have every m> d, 


SmiU^)p > 2-^‘=‘/3(d,TO)2-“'" > ( 2 “-l)-i 2 -^^ 2 -“"‘ 
and therefore, 

Sm(C/“)p > ( 2 “ - l)-i 2 -^^(d- 


m 
d-l)' 


( 2 . 12 ) 


(2.13) 


Proof. The inequality (12.131) follows from (12.121) and the inequality (^(0^) > (j^)^ ^prove (I2.12p . Let 
Mi be the cubic cardinal B-spline with support [0,4] and knots at the points 0,1, 2, 3,4. Since the support of 
functions M4(2^+^-) for k € Z_|_ is the interval [ 0 , 2 ^], we can extend these functions to an 1 -periodic function 
on the whole R. Denote this periodic extension by pk- Define the univariate nonnegative functions gk.s and 
on T by 


9k. 


,s{x) ■■= fj,k{x-4:s), fc e Z+, s = 0, ...,2'" - 1; g^{x) := gk,s{x). 


(2.14) 


s =0 


From the identity M 4 ^^(x) = M 2 (x) — 2 M 2 (x — l)+M 2 (a: —2), X € R, it follows that |M 4 ^^(x)| < 2 M 2(2 ^x), xG 
R, and consequently. 


\9l^],(.x)\ < 2'^'"+^pk,s{x), xeT. 

One can also verify that 

suppgfe,s = Ik,s =: [2“'"s,2“''(s + 1)], int (~l int 4,s' = 0, s ^ s', 
and by the equation Mi{x) dx = 1, 


/ gk,six)dx = 2 
Jo 


-k-2 


Let 


9k,s ■ 


,s := kGZl, sGl\k). 


1=1 


Put for u C [d], := T=p—-gyr- Then we have for k G s G I'^{k) and u C [d], 

1 Lj£u '^j 

\Dl9k,six)\ = 2^^''^2‘^^'^'^^^Y[Tk,i.sjixj)Y[9k,i,siixj),xGT‘^. 

jeu j^u 

Take m, n S N with n > m and define the function 

n 

fm,n-.= 2-5'^^] 2-“' y] y] gk,s. 

l=m \k\i=l sGl'^ik) 

Observe that 

n 

fm,n= 2-5^'^] 2-“' y] gk, 
l—m \k\i—l 

where 


(2.15) 

(2.16) 

(2.17) 


(2.18) 


(2.19) 


9k ■■= ^ 

1=1 


( 2 . 20 ) 
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We will show that 


jeu 


( 2 . 21 ) 


Let us prove this inequality for u = [d] and h G , the general case of u can be proven in a similar way with a 
slight modification. For h G K and fc G N, let h^^'> G [0,2“^) is the number defined by h = + s2“^ for some 

s G Z. For h & and k G put 

hW := 

1=1 

Since gk is 2“*^i-periodic in variable Xj, we have 

^h^'^\9k,x) = All^}{gk,x), xGT'^, hGW^. 

By using the formula 

^hif,x) = [ {x + y)[h-^M2{h-^y)]dy 

Jr 

for twice-differentiable function / on K (see, e.g., [5J p.45]), and (12.151) we get 



and therefore, 

^h^'^\9k,x) = ^hmi9k,8,x) 


s&I'^ik) 


= / ^[d]9kAy)'[[ilT-f^) \yj-Xj))dy. 

s6/d(fe)j=l j=l 


Hence, 


|A"’‘"k5fe,a;)l < E / 5A\>^'AkAy)t[ihf^ft[ihf)-^M2iihf'>)-\y,-x,))dy 

se/‘i(fc) 1=1 1=1 

SG/-I(fc) 1=1 1=1 1=1 

< 25'^2l'=l^n^“ / E ‘PkAy)l[(hf)-^M2{(hf)-\yj-Xj))dy. 


' se/‘*(fc) 


By the inequalities (12.181) and 


E E E V’k,s{y) < 1, y G M'', 

l—m |fc|i—/ sG/‘^(fc) 
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we derive that 


l—m \k\i—l 



Xj))dy 


d p d d 

= !!'■?/, s n*?- 

j=i j=i j=i 

The inequality (12.211) is proven. This means that fm,n G U^. From (12.201) . (12.161) . (12.171) we have that if m is 
given then for arbitrary n > m, 


n ” // _ 1 \ 

WUnWp > WUnWi = 2-5'^ ^2-“' ^ \\g,\U = J- (2-22) 

l—m |fc|i=/ l—m ^ ' 

On the other hand, one can verify that fm,,niO = 0) C G G(m) which yields that fm,n) = 0 for arbitrary 

$m, and consequently, by ( 12 . 221 ) for arbitrary n > m, 


Sm{UZ,)p > \\fm,u-SU<i>Jm,u)\\p = \\U.n\\p > 2“"" E 2““'^ _ M 

l=m ^ ' 

This means that for t = 2““, 

Therefore, applying (12.111) gives 

s™(f7“)p > 2-7''2-“’"tF™,d-iW = 2-7''2-“™/3(a,d,m) 

> 2 "'^‘' 2 ““’"( 2 “ - 1 )"'^ ( 2 “ - l)"^"^ 

= ( 2 “ - l)-i 2 "'^‘* 2 "“™ 

which proves (I2.12F □ 

Corollary 2.1 Let 0 < p < oo and 0 < a < 2. Then we have for every m > d, 
2-7''/3(d,m)2-“’" < Sm{U^)p < [2(p+l)i/P]-‘^/3(d,m)2-“™. 
Moreover, if in addition, m > 2{d— 1), 

(2a_i)-l2-7d2-am|^^^ J < 5^(17“ )p < a%d) 2-^ . 
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Put 




Theorem 2.5 Let 1 < p < oo, 0 < a < 2 and 1 < v < d. Then we have every m> v, 

s'LiU^np > 7'(^,m)2— > [2^(2“ - I)]"! (^2“ - ' 2 — J. 


(2.23) 


Proof. We first consider the case v = d. We take the univariate functions gk,s: Ok as in (I2.14L and for every 
u C [d], define the functions 

9k,s ■= ndfej.sj, A: G Z(J., s&I'^{k), g'k'-= 

jeu iGu 

If m, n G N with n > m, and u C [d], we define the |u|-variate function 


:= 2-5|“I ^ 2 -“' ^ ^ 

l—m / sGl^{k) 

and the d-variate function 

'pTn,n ■= fm,n- 

uC[d] 

From the proof of Theorem 12.41 we can see that 
f>m,n G U^, 4>m,ni£,) = 0, ^ G G'^ {m). 


(2.24) 


Observe that 


:= ^ ^ 2 -“' ^ gl 

uC[d] l—m |fc|i=Z 

From (12.201) . ()2.16|) . (12.171) we can also verify that if m is given then for arbitrary n > m, 

n 

UmAv > Um,nh = E 2-5H ^ 2““' E 

uG[d] l—m |fc|i=Z 


= 2“^'“' E 2 


l=m ^ s=l ^ l=m 


— al 


1-1 
S — 1 


By (12.241) we have = 0 for arbitrary 4)^, and consequently, by (12.221) for arbitrary n > m, 

d 


stiUAp > UmAp > EQ2'"^E2-“'P_^V 

s—1 ^ ^ l—m ^ ' 


(2.25) 


This means that for t = 2 




S =1 




l—m 


S — 1 


s=l 


(2.26) 
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Therefore, applying (12.31) gives 


> 2 "“™ J2 J2 ( 7 ) 

s=l i =0 ^ ' 

E C+ 1 ) E ( 7 ) -1)' 

s^O V / /^o ^ ^ 

W (2“ - 1)' E C + 1 ) (2“ - 1)-^^+'^ 

1^0 ^ ^ s^O ^ ^ 

> 2"“™ (2“ - l)^“i E ~ (2“ - l)-("+i) 

= [2^(2“ - l)]-i 2"“™ J (2“ - l)-^-! E (^ ~ ^) [2^(2“ - 1)]“" 


= 2 “ 


= 2 ’ 


= [2^(2“ - 1)]-^ (2“ - 1)'^-! [1 + 2^(2“ - 1)-2-“™ 


m 

d-1 


which proves (12.2311 for the case v = d. 

To process the case u < d we take the function (j)m,n G C/S" by 

M V p 

^m,n / j Jm^n 

u<Z[v] 

Similarly to (12.251) we have 


> ikl^!„iii = E(j2-"^E2 


■\—al 


i - 1 
S — 1 

s=l l—m 

Hence, by replacing d by z/ in the estimations (12.261) and (12.2711 we prove the case < d. □ 
Corollary 2.2 Let 0 < p < 00 , 0 < a < 2 and 1 < v < d. Then we have every m> v, 


1=0 


2 ”“™E ( 7 ) < 2 -“™^ 


/=o 


i/p]-i 


13(1,m). 


Moreover, if in addition, m > 2(v — 1), 


2Y“-1)]-M2“-^)" ' 2 —f <sS(C/Sa<«Y2 


2.5 Cubature formulas 


We are interested in cubature formulas on Smolyak grids for approximately computing of the integral 
Hf)~f f(x)dx. 


(2.27) 
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If / G we use the cubature formula on grids G{m) given by 

LiAmJ) = Y. 

^eG(m) 

where Am = integration weights. The quantity of optimal cubature Intm(Td) on Smolyak grids 

G(m) is defined by 

IntmiFd) := inf sup |/(/) -/)|. 

Am feFd 


For a family = {¥’€}^gG(m) functions on the linear sampling algorithm •) generates the 

cubature formula /) on Smolyak grid G{m) by 

im{AmJ) = Y 

where the integration weights A^ = (\)^gG(m) given by = fj^(p^(x) dx. Hence, we have |/(/) — 
im{Am,f)\ < 11/- 5r„($m,/)||i and consequently, 

Intm(Fd) < °Sm{Fd)i- ( 2 . 28 ) 


Theorem 2.6 Let 0 < a < 2. Then we have for every every m > d, 
2 -^‘^^(d,m)2-“’"^(d,m) < InC(t/“) < 2-2'^/3(d, m) 2-“"*. 
Moreover, if in addition, m > 2{d — 1), 

(2“-l)-i2-7^2-“'" J < ikmiU^) < a°{d)2-^‘^2-^^ 

where 


{ 1, a > 1, 

d, 0=1, 

(2 “-l)-rf, a<l. 


(2.29) 


(2.30) 


Proof. The upper bounds in (I2.29l) " (l2.30p follow from (12.281) and Theorem 12.21 To prove the lower bounds, 
we take the function fm,n G as in (12.191) with the property fm,n{f,) = 0, / G G{m). Notice that fm,n is 
a nonnegative function. Hence, we have A^($m,/m,n) = 0 for arbitrary and consequently, by (12.221) for 
arbitrary n > m, 

InU(17“) > |/^.„-A^(c!>^,/„.„)! = Wfm.nWi > 

Comparing with (12.121) . we can see that Int,^(17^) can be estimated from below as in the proof of Theorem 
(12.41) for Sm{U^). This proves the lower bounds in (I2.29p - (I2.30I) . □ 

If / G U^'^, we use the cubature formula on grids G‘'{m) given by 

c(a:;,/) = E 

(m) 

where A'f,^ = are integration weights. The quantity of optimal cubature IntmiFd) on Smolyak grids 

G'^{m) is introduced by 

InC(Fd) := inf sup |/(/) -/^^(A^,/)|. 

A- 
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Put 




;=o 

In a similar way to Theorem 12.61 we obtain 

Theorem 2.7 Let 0 < a < 2 and 1 < ly < d. Then we have for every every m > v, 

< InC([/“.-) < 2““™ ^ 2-2'/3(|, m). 

1=0 ^ ' 1=0 ^ ' 


Moreover, if in addition, m > 2(z^ — 1), 


[ 2 ^( 2 “ - 1 )]-^ 2 “ - 


127 

m 


ty-i 


V- 1 


< InC(t/“-'^) < a{iy)T 


u-ir 


where 


(5/4)^ 

a > 1, 

1^(51 AY, 

a = 1, 

[1 +1/4(2“-l)]^ 

a < 1. 


3 Sampling recovery based on B-spine quasi-interpolation 


In this section, we extend the results on upper bounds of s'f^(U^^)p and in Subsection l2.dl to arbitrary 

smoothness a > 0. To this end we construct B-spline quasi-interpolation representations for continuous functions 
on find an explicit formula for the coefficients functionals and estimate the quasi-norm of Lpiff’^) for the 
component functions in these representations for functions in i7^(T^). 

3.1 B-spline quasi-interpolation representations 

In order to construct B-spline quasi-interpolation representations for continuous functions on we prelimi¬ 
narily introduce quasi-interpolation operators for functions on For a given natural number £, denote by Mi 
the cardinal B-spline of order I with support [0,i] and knots at the points 0,1, We fixed r S N and take 
the cardinal B-spline M := M 2 r of even order 2r. Let A = {A(s)}|j|<^ be a given finite even sequence, i.e., 
X{—j) = X{j) for some yi > r — 1. We define the linear operator Q for functions / on R by 

Q{f,x) ■■= '^Hf,s)M{x-s), (3.1) 

sGZ 

where 

Mf,s):= ^ A(j)/(s - j-f r). (3.2) 

bl<M 

The operator Q is local and bounded in C'(R) (see [H p. 100-109]). An operator Q of the form (I3.1I1 - (|3.2I1 is 
called a quasi-interpolation operator in C'(R) if it reproduces V 2 r-i, i-e., Q{f) = / for every / € V 2 r-i, where 
Vi denotes the set of d-variate polynomials of degree at most Z — 1 in each variable. 

If Q is a quasi-interpolation operator of the form (ICT-lT^. for h > 0 and a function / on R, we define 
the operator Q(-; h) by Q{f; h) := ah o Q o ai/h(f), where ah{f, x) = f(x/h). Let Q be a quasi-interpolation 
operator of the form in C(R). If fc € Z+, we introduce the operator Qk by 

Qk{f,x):= Q(/,cr;hW), x e R, := {2r)-^2->^. 
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We define the integer translated dilation Mj-^s of M by 
Mk,s{x) '■= M{2r2^x — s), k € s G Z. 

Then we have for k G 

Qk{f){x) = '^ak,s{f)Mk,s{x), VxGR, 

sSZ 

where the coefficient functional ak,s is defined by 


afe,s(/) := A(/,s;/r('=)) = X{j) f - j + r)). (3.3) 

bl<M 

Notice that Qk{f) can be written in the form: 

Qk{f){x) = + r))Lfe(x - s), Vcc e R, (3.4) 

sez 

where the function is defined by 

Lk:== MJWk,,- (3.5) 

bl</^ 

From (13.41) and (13.51) we get for a function / on R, 

IIQfc(/)llc(R) < I|Aa||c(r)II/IIc(r) < I|A||||/||c(r), (3.6) 

where 

La(x) := = E E Mj)M(x-j-s), IIAll = |A(j)|. (3.7) 

sez|j|</j b'l</i 


For k G Z^, let the mixed operator Qk be defined by 

d 

Qk-=X{Qk,, (3.8) 

i=l 

where the univariate operator Qk^ is applied to the univariate function / by considering / as a function of 
variable Xi with the other variables held fixed. We define the d-variable B-spline Mk^s by 


d 

Mk,s{x) := Y\_Mk,,si{xi), kGZ'l, sGZ'^, (3.9) 

i^l 

where Z'^ := {s G Z’^ : Si > 0, i G [d]}. Then we have 
Qkifyx) = ^ ak,sif)Mk,six), VxGR'^, 

seZ‘1 


where Mk^s is the mixed B-spline defined in (13.9L and 


akAf) = (Y[<^kj,s,']if), 

\=i / 


(3.10) 


and the univariate coefficient functional aki,si is applied to the univariate function / by considering / as a 
function of variable x, with the other variables held fixed. 
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Since M{2r2^x) = 0 for every k G Z+ and x ^ (0,1), we can extend the univariate B-spline M(2r2^-) to 
an 1-periodic function on the whole K. Denote this periodic extension by Nk and define 

Nk,s{x) ■■= Nk{x-s), fc G Z+, sGl{k), 

where I(k) := {0,1,2r2^ — !}■ We define the d-variable B-spline Nk^s by 

d 

NkAx):= sGA{k), 

where I'^{k) := 0^=1 1{ki)- Then we have for functions / on 

Qkif,x) = ^ ak,sif)Nk,six), VxGT"^. (3.11) 

sGl^{k) 

Since the function L\ defined in (1211) is 1-periodic, from (13.61) it follows that for a function / on T, 

IIQfc(/)llc(T) < ||Ta||c(t)||/||c(t) < I|A||||/||c(t), (3-12) 


For k G Z'l, we write fc —>■ oo if > oo for i G [d]). In a way similar to the proof of [141 Lemma 2.2] one 
can show that for every / G C'(T‘^), 

Wf-QkimCiT^) < C T. 

uC[d], u^0 


and, consequently, 

11/ “ Qfc(/)llc(T‘i) 0, fc —>■ oo. 


(3.13) 


For convenience we define the univariate operator Q-i by putting Q-i{f) = 0 for all / on I. Let the 
operators qk be defined in the manner of the definition (13.81) by 


i=l 


Qk ■■= W[Qki-Qki-l\ fc G Z+. 


(3.14) 


From the equation Qk = J2k'<k 9^' p.l3l) it is easy to see that a continuous function / has the 
decomposition / = Qkif) with the convergence in the norm of C(T‘^). From the refinement equation 

for the B-spline M, in the univariate case, we can represent the component functions qk{f) as 


'?/=(/) = X! ^k,sif)Nk,s, 

sGl'^ik) 


(3.15) 


where Ck,s are certain coefficient functionals of /. In the multivariate case, the representation (13.151) holds true 
with the Ck,s which are defined in the manner of the definition (13.101) by 

CkAf) = 

\=i / 

See [13] for details. Thus, we have proven the following periodic B-spline quasi-interpolation representation for 
continuous functions on T"^. 


Lemma 3.1 Every continuous function f on is represented as B-spline series 

/ = E 9fe(/) = E E ^kAf)Nk,s, 


(3.16) 


converging in the norm of where the coefficient functionals CkAf) explicitly constructed as linear 

combinations of at most mo of function values of f for some mg G N which is independent of k,s and f. 
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3.2 A formula for the coefficients in B-spline quasi-interpolation representations 

In this subsection, we find a explicit formula for the coefficients Ck,s{f) and hence, estimate the quasi-norm of 
|| 9 fc(/)||(? of the component functions in the periodic B-spline quasi-interpolation representations (lfi.1611 . 

If h G we define the shift operator for functions / on by Th{f ) := f{- + h). Recall that a d-variate 
Laurent polynomial is call a function P of the form 

P{z)= (3.17) 

seA 

where Al is a finite subset in and z® := n,ti . A d-variate Laurent polynomial P as (13.171) generates the 

rpi 

operator Tj^ by 

^^r(/)= (3.18) 

sSA 

fPl \P(z)] 

Sometimes we also write ^ l^otice that any operation over polynomials generates a corresponding 

fPl 

operation over operators . Thus, in particular, we have 

= aiTf^'(/) + a 2 Tf^'(/), o (/)• 

By definitions we have 

A' = A := f[izj - 1)', A,. := Y[(.Zj - !)'• 

j=l jGu 

We say that a d-variate polynomial is a tensor product polynomial if it is of the form P{z) = 
where Pj{zj) are univariate polynomial in variable Zj. 

Lemma 3.2 Let P he a tensor product Laurent polynomial, h G with hj ^ 0, and I G N. 

Tjf'\g) = 0 for every polynomial g G A-i Then P has a factor Di and consequently, 

= T’r'oAi, p{z)=DiP*{z), 
where P* is a tensor product Laurent polynomial. 


rii^i 

Assume that 


Proof. By the tensor product argument it is enough to prove the lemma for the case d = 1. We prove this case 
by induction on 1. Let P{z) = Yl'^=-m^sZ^ for some m,n G Z+. Consider first the case I = 1. Assume that 

Plf^g) — 0 for every constant functions g. Then replacing by go = 1 in (13.181) we get T^^^(go) = J2^=-m ~ 

0. By Bezout’s theorem P has a factor (z — 1). This proves the lemma for 1 = 1. Assume it is true for 1 — 1 and 

rpl 

A ( 5 ) “0 for every polynomial g of degree at most 1 — 1. By the induction assumption we have 

° A'-i, P(z) := (z - l)'-'Pi(z). (3.19) 


We take a proper polynomial gi of degree 1 — 1 (with the nonzero eldest coefficient). Hence ipi = AJ^ ^{gi) = a 
where a is a nonzero constant. Similarly to the case 1 = 1, from the equations 0 = Tjf’\gi) = we 

conclude that Pi has a factor (z — 1). Hence, by (13.191) we can see that P has a factor (z — 1)^ The lemma is 
proved. □ 


Let us return to the definition of quasi-interpolation operator Q of the form induced by the sequence 
A as in (13.21) which can be uniquely characterized by the univariate symmetric Laurent polynomial 


Pa(z) := z’' E A(s)zL 
|s|</i 
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Let the d-variate symmetric tensor product Laurent polynomial Pa be given by 


Pa(z) := ^ A(sj)zL 

j=l 

For the periodic quasi-interpolation operator Qk(f) = o-k,sif)Nk,s given as in (13.1 111 , from (13.311 we 

get 

akAf) = T‘f^)’(/)(sh«). (3.20) 

Let us find an explicit formula for the univariate operator qk{f)- We have for k > 0, 

Qkif) = E 

se/(fc) 

= E Ti%iimsA'^^))Nk,2s + E Ti%iifms+i)AA)Nk,2s+i. 

From (13.2011 and the refinement equation for M, we deduce that 

Qk-i(f) = E 

sG/(fe-l) 

j=o A sei(k-i) 




+ 2-^’'+^E(2, + i) E 

j=o \ J ^ y sei{k-i) 

=■■ QT-iiD + Qk-Mf)- 

By the identities Nk^ 2 r 2 '‘+m = Nk,m and f {A^">){2r2^ -\-m) = /{A^^rn) for fc G and m G Z, 

we have 

Qr-T(/) = 2-2’'+'E(2") ^ T™(/)(2 (s-j>«))7V,,2s 

j=o A y sGi+/(fc-i) 

= 2-"’'+'E(2") E r™(/)(2(s-J>('=)))iVG2. 

i=o A y sei(k-i) 

sG/(fe-l) 

where 

P'_(z) := 2-^^+^PAz^)jZ ® (3-21) 

j=o A 

In a similar way we obtain 

QtAu) = E 7l{^r'(/)((2s + i)/i('=^)iVM.+i, 

se/(fc-i) 

where 

PUAz) ■■= 2-2’'+1Pa(z2)2 ( 2/+ (3.22) 
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We define 


Peven Pa Pevem Podd ■— Pa -Podd (3.23) 

Then from the definition qkif) = Qkif) — Qk-iif) we receive the following representation for qkif), 

'?o(/) = E (3.24) 

se/(o) 

and for /c > 0 , 

%(/) = 9r"(/)+9fc""(/) (3.25) 

with 

sGl(k-l) 

s 6 /(fc-l) 


9r"(/) = 

9^ (/) = 


From the definitions of Qk and qk it follows that 
r^f;E( 5 )( 2 s/i('=^) =0 and (g)(( 2 s + l)/rW) = 0 for every g€Vr. 

Hence, by Lemma 13.21 we prove the following lemma for the univariate operators qk- 


Lemma 3.3 FFe have 

Peven{z) = i42r (z)Peven(^) 

Podd(z) = D2r{z)P:^d{z). 

where Peven > Podd ® Symmetric Laurent polynomial, 
we have for fc > 0, 


(3.26) 

Therefore, in the representation (13.241) - (13.251) of qk[f ), 


9 r"(/) = E ° {msh^^))Nk, 2 s, 

sG/(fc-l) 

= E ° A2J.,(/)((2s + l)/rW)7Vfc,2.+i. 

sG/(fe-l) 

Equivalently, in the representation ()3.15p of qk{f), we have for s G 7(0) 


coAf) = Tl%]\f)isAA, 

and for k > 0 and s G 7(fc), 

s even 
s odd. 


Tl^ir^oAl\Af){sAA, 
AJf-A Al%,if){sAA, 


Proof Consider the representation (13.151) for qk{f) and d = 1. If g is arbitrary polynomial of degree at most 
2r — 1, then since Qk reproduces g we have qk{g) = 0 and consequently, CkAA = 0 for fc > 0. The equations 
(I3.24|) - ()3.25l) give an explicit formula for the coefficient Ck,s{g) as and Tj^^fA\g){(2s + l)A'^'>). 

Hence, by Lemma 151^ we get p.26l) . □ 
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(3.27) 


Theorem 3.1 In the representation H3.15I) of qk{f), we have for every k G Z+(u) and s G I‘^{k), 

ckAf) = 

where 

PkAA = n ) n (^.-) n ^2. ), 

j^u jGu jeu 

f Pe^renA)^ « even, 

Pk,,sMl) = { 

1 -Podd(2j). s odd- 


(3.28) 


(3.29) 


Proof Indeed, from the definition of Ck,s{f) and Lemma 1531 we have for every k G Z^|.(m) and s G A{k), 

CkAf) = = Tl^A\f){sAA, 

S=i ^ 

where 

(PkAj), kj=0, 

Pk^,sAzj) = lPf^^Azj)D 2 rizj), kj>0, seven 

[pAdAj)D 2 r{zj), kj > 0, s odd. 

□ 


For a Laurent polynomial P as in (13.171) we introduce the norm 

ll^lh= Elc«l- 

sgA 

Notice that ||La||c(t) < ||-Pa||- 

For given A, 1 < p < oo and 0 < a < 2r, we define 

a = a{p,a,A,r):= ( 2 r)-“[r(p + 1 )]"^/^ max {||Pj; 3 „||, ||P„*dd||}, 

(3.30) 

b = b{A) := ||La||c(t). 

Theorem 3.2 Let 1 < p < oo, and 0 < a < 2r. Let f G Then f can be represented by the series 

(I3.16|) converging in the norm of C{T‘^). Moreover, we have for every k G 'Z'^{u), 

\\qk{f)\\p < al“l 6 "-l“l 2 -“l'=lM/|//»(.). (3.31) 

Proof. The first part of the lemma on representation and convergence is in Lemma l3.II Let us first prove p.31l) 
for p = oo. For convenience, we temporarily use the notation a = Qp. In this case, (13.311) is as 

||gfc(/)||oo < sup \ck,sif)\ < al(^l 2"“l*li l/l/f. („). (3.32) 

sei'^(k) 

By definition we have 



Hence, by (I3.12|) we derive 

lkfc(/)||oo < 6'^-'“' \\g\A, (3.33) 
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Similarly to the proof of Theorem I2.1l we can show that for every k G Z^(u), 


II5II00 = 

(l[<lk\f) 

< sup 

fn 


\ J 

00 se/‘'(fc) 

V / 


Observe that for every Laurent polynomial and every h G M, 

iiTr(/)iioo < iiPiiii/iioo. 

Setting 

jeu 


(3.34) 


(3.35) 


we get 



T, 


K,s 

hm 


h(^) 


(/,shW) 


Hence, by (|3.27I) ~ (I3.29I) and (13.351) 


sup 

fn 

= sup 

rr^lPLj 

A',r.)(/,shW))j 

sGl^{k) 

V / 

sGl^{k) 


L -1 


< sup A^“)(/,s/i(''))) 

sei‘‘(k) 

< [max{||P;_||,||P;ddll}]'“' ni^f 

jeu 

This together with ()3.33|) and (13.3411 proves (13.3211 . 

Let us prove (13.3111 for the case 1 < p < oo. We have for every k G Z'^{u), 

p 

^k,s{f )^^k,s (^) 

se/‘*(fc) 

> ; Mk,s{x)' 

Jfd 


lkfc(/)llp = / CkAfWkA^) 


dx 


< sup |cfc.s(/)|^ 
se/'^(fe) 


Mk,s{x) 

sGl‘‘(k) 


dx 


= sup \ck,s{f)\^\I'^{k)\ [ \Mk,o{x)\Pdx 


= sup |cfc.s(/)r|I'^(/c)| TT / \M{2’^^Xj)\P dx^ 

s&Id^ik) 

= sup |cM(/)n/yA:)|(2r)-H2-l'=li f / |M 2 ,(t)rdt 

sG/‘^(fe) \JK 

By employing (13.3211 . Young’s inequality 

\\M2r\\Lp(S.) = \\M 2 r -2 * M2\\lp(R) < ||-W2r-2 || Li (R) ||-^2 || Lp(R), 


l“l 
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the equations ||M 2 r- 2 ||Li(B) = 1 and ||M 2 ||^ 


= 2{p + 1) we complete the estimation as follows 


hkifWp < 


^\u\ ^d-\u\ 2 -a|fe|l 


P 






= Kp+ir'"' 




2-pa\h\i 


I/I 


P 




2-pa\k\i 


I/I 


P 


□ 


Theorems on B-spline quasi-interpolation representations with discrete equivalent quasi-norm in terms of 
coefficient functionals have been proved in m-m for non-periodic various Besov spaces. 


3.3 Sampling recovery and cubature 

For m € Z+, we define the operator by 

^™(/):= E = E E ^m(/)^M- 

: I supp(fc) I <(/, \k\i<m : | supp(fc) | <1^, |fc|i<m sGl^{k) 

For functions / on having at most i/ of active variables, defines the linear sampling algorithm on the 
Smolyak grid G^{m) 

K.{f) = = E 


where n := |G‘"(m)|, := {V'c}{eG''(m) and for ^ = 2 ^s, ip^ are explicitly constructed as linear combinations 

of at most at most N B-splines Nkj for some TV S N which is independent of k, s, m and /. 


With where a, b are as in (Id.dOll we put 


(5(z/, a, b, m) := E! 

1=0 

and 


b’' ^ 


c{iy,a,b):= |2“ 


2|-| sgn(a-l)| ^ 


[a + bY, 

a > 1, 

iy[a + bY, 

0 = 1, 

[«/(2«-l) + 6]-, 

a < 1. 


Theorem 3.3 Let 1 < p < oo, 0 < a < 2r and 1 < v < d. Then we have for every m > v, 


< sup \\f - Rfr,if)\\p < 5{v,a,b,m)2- 
f&uSd’' 


(3.36) 


1 / n-arn, 


< exp(2“ - 1) [a/(2“ - 1) -b b^ 2 
Moreover, if in addition, m > 2{y — 1), 


s'LiU^lp < sup \\f - R‘fr,if)\\p < c{v,a,b)2- 

fauSd'' 


i/- 1 


(3.37) 
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Proof. We prove the theorem for the case u = d which yields the case v < dhy the same argument as in the 
proof of Theorem 12.21 Let / € Put t = 2““. From Theorem KL2I and and (12.41) it follows that for every 
m > d — 1, 

Wf-Riimp < E E 

^\u\ i^d-\u\ 2-a|fc|i 

|u|<rf |fc|i>m 



d 


1^0 

d 




; |A;|i>m 

j —ct(m+j) 

l-l 


/=o ^ ^ i=i 

= 2-“(™+l) E (f) 
1=0 ^ ^ 


_ 2 °‘’^S{a,d,a,b,m). 

The first inequality in (13.361) is proven. 

Hence, applying Lemma [2.31 for t = 2““ gives 


Wf-Riimp < exp(2“-l)2-“-E(?)“'^"”'[2“-l]“'^'-' 

1^0 ^ ^ 

< exp(2“ - 1) [a/(2“ - 1) + bf 2““'" 
which proves the second inequality in (13.361) . 

We now prove (13.371) . If m > 2{d — 1), by (13.381) and Lemmafor t = 2““ we get 

\\f-Rt{f)\\p < 2-“(-+i)E(;)a'&"“'(;ri)^'-i(2-“) 


(3.38) 


< 2 


-a(m+l) 


m 

d-l 


d 

E 

1=0 


a' 6“‘-'6,_i(2-“) 


= c(a,d,a,&) 2 -“-(^^^ J. 

The proof is complete. □ 

In a similar way to Theorem 12.61 from Theorem 13.31 we obtain 

Theorem 3.4 Let 0 < a < 2r and 1 < v < d. Let a = (2r)“““^ max {||Pgyg„||, ||Poddll} ^ = l|iA||c(T) (cf- 
(I3.30|) ). Then we have for every m > u, 

InC(17“’‘^) < 6{v,a,b,m)2-°^'^ < exp(2“ - 1) [a/(2“ - I) + 6]“^ 2-“"^ 771 “^-^ 

Moreover, if in addition, m > 2{v — 1), 

InC(C/“’-) < c(77,a,6)2—J. 
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4 Examples 


The upper bounds for s'^{U^'')p and Int5:^(?7^’^) obtained in Theorem 13.31 and Theorem 13.41 depend on the 
parameters p, a, d, v, r, a, b. In applications, if the parameters p, a, d, v are a priori known, these upper bounds 
are controlled by parameters r, a, b which are determined by the choice of a univariate quasi-interpolation 
operator Q of the form (EH). The operator Q is induced by the sequence A as in (III.2I1 which can be uniquely 
characterized by the univariate symmetric Laurent polynomial Pa- Moreover, the parameters a,b defined in 
dSHOl) contain max {||P;^g„||, ||P*ddll} and ||La|1c(t) which are desirable to be minimum by the choice of Pa- 
In this subsection, we give some examples of the univariate symmetric Laurent polynomial Pa characterizing 
quasi-interpolation operator Q of the form EH) with approximate estimates of ||Pevenll! l|.Po*ddll I|.^a|Ic(t) 
or ||Pa|| with ||La||c(t) < I|.Pa||. For a given Pa, the Laurent polynomials Peven ^’^cl P*^^ can be computed 
from and (ESH- 

4.1 Piece-wise linear quasi interpolation 

Let us consider the case r = 1 when M{x) = (1 — |a; — 1|)+ is the piece-wise linear cardinal B-spine with knot 
at 0,1, 2. Let A = {A(s)}j=o (a* = 0) be a given by A(0) = 1. If Nk is the periodic extension of M(2*^+^.), then 

Nk,s{x) ■-= Nk{x-s), fc e Z+, se/(fc), 

where I{k) := {0,1,..., 2*+^ — 1}. Consider the related periodic quasi-interpolation operator for functions / on 
T and k G 

Qk{f,x)= Y. /(2-("+')(s + l))iVfe,,(x) 

s6/(fc) 

We have 

Pa{z) = z, 

Peven{z) = — — (z—1) , , 

Poddiz) = Poddiz) = 0 , 

and 

iifaII = 1 , iiFe;e„ii = \\p:dd\\ = «. 

Hence, 

1 

doif) = ^T™(/)(2-'s)iVo.s = /(o), 

s^O 

and for /c > 0, 

dkif) = C“(/) = 

With these formulas for qk{f) after redefining Nk^ 2 s as ipk,s, the quasi-interpolation representation (13.161) be¬ 
comes the Faber series. 


4.2 Cubic B-spline quasi-interpolation 

For r = 2, we can take 

Pa{z) = ^(-z + 8-z"i), = + -^z- 

O ODD 
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Then, we have 


Peveniz) = {z - 1 )* P*^^n{z), P*^en{z) = ^Z ^ (^z"^ + Az^ + 8 z^ + Az + , 

-Podd( 2 :) := {z — l)^Po*(i(j(z), T’o*dd(^) ■= Y2 + 4 z + , 

and 

II^aII = y, l|f’;ve„ll = I ll^o*ddll = i- 

4.3 Quintic B-spline quasi-interpolation representation 

For r = 3,we can take 

PAiz) := Y4400 
Then, we have 

PevenW = {z - if P:,,,{z) , Podd(z) = (z-lfP^VW, 


Peve„(^) = [139760+ 97002(z + z-i)+42508(^2+ z-2) + 11462(z3 + z-3) 

+ 2328(z^ + z-f + 458(z'^ + z"*^) + 36(z® + z“®) + 6(z^ + z"’')], 

^odd(^) = 77J7^ [164910 + 97002(z + z-f + 36632{z^ + z-^) + llA62(z^ + z-^) 

4oUoUU 

+ 2776(z^ + z"^) + 458(z'^ + z"*^) + 88(z® + z“®) + 6(z^ + z"’’) + (z® + z"®)]. 


and 

II^aII 


37904 

14400 


2.63, 


IP* 


447360 

460800 


0.97, 


IIP 


ddl 


467160 

460800 


1 . 00 . 
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